Scalar Mesons a (1450) and a(600) from Lattice QCD 
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We study the ao and a mesons with the overlap fermion in the chiral regime with the pion mass 
as low as 182 MeV in the quenched approximation. After the rj'rr ghost states are separated, we find 
the do mass with qq interpolation field to be almost independent of the quark mass in the region 
below the strange quark mass. The chirally extrapolated results are consistent with ao(1450) being 
the ud meson and A"q(1430) being the us meson with calculated masses at 1.42 ± 0.13 GeV and 
1.41 ±0.12 GeV respectively. We also calculate the scalar mesonium with a tetraquark interpolation 
field. In addition to the two pion scattering states, we find a state at ~ 550 MeV. Through the 
study of volume dependence, we confirm that this state is a one-particle state, in contrast to the 
two-pion scattering states. This suggests that the observed state is a tetraquark mesonium which 
is quite possibly the <r(600) meson. 
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I. INTRODUCTION 

Unlike pseudo-scalar, vector, and tensor mesons, 
the scalar mesons are not well known in terms of 
their SU(3) classification, the particle content of 
their composition, or their spectroscopy. Part of 
the problem is that there are too many experimen- 
tal candidates for the qq nonet. Fig. Q] shows the 
current experimentally known scalar mesons whose 
number more than doubles that of a nonet. One vi- 
able solution is that low- lying scalars, such as the 
er(600), ao(980) and /o(980) are tetraquark mesoni- 
ums whose classification and spectroscopy have been 
studied in the MIT bag model [l[ and the poten- 
tial model 0. Another suggestion is that ao(980) 
and /o(980) are KK molecular states [3|. Other 
candidates for tetraquark mesoniums include vec- 
tor mesons pairs produced in 77 reactions [1] and 
hadronic productions Q and the recently discovered 
charmed narrow resonances @ . 

Under the supposition that ao(980) and /o(980) 
are tetraquark mesoniums on account of the fact 
that they are favored by spectroscopy studies [l], 0] , 
small two-photon decay widths Q , and the pattern 
of <f> and J/*f? decays |8[ , the question remains: where 
is the isovector scalar qq state? From Fig. [TJ we 
see that one candidate is do (1450). However, in the 
conventional wisdom of the quark model, its mass is 
too high. Not only is it higher than 02 (1320) and 
ai(1230), in contrast to the spin-orbit splitting pat- 
tern in charmonium; it is even slightly higher than 
i^o(1430) which contains a strange quark and is be- 
lieved to be the su or sd meson in practically all 
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FIG. 1: Spectrum of scalar mesons together with 7r, p, a\ 
and <J2 mesons. 



the models [9J. According to the quark counting 
rule, mesons and baryons made up of strange quarks 
are expected to lie higher than their counterparts 
with u/d quarks. Notwithstanding the success of the 
quark potential model in describing charm and bot- 
tom hadrons, its applicability to light hadrons with 
SU (6) symmetry has been questioned, since chiral 
symmetry plays an essential role [13, El in light 
hadron dynamics. Might it be that the scalar qq 
meson is yet another challenge to the SU(Q>) quark 
model's delineation of light hadrons? 

Lattice QCD is perhaps the most desirable tool 
to adjudicate the theoretical controversy surround- 
ing the issue and to reveal the nature of the scalar 
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mesons. In fact, there have been several calcula- 
tions to study the ao meson with the tptp interpola- 
tion field in the quenched approximation fl2l . fl3l [l4| 
and with dynamical fermions [HI, EH, E3, EH, EH • In 
calculations with relatively small quark masses, it is 
found that the ao mass does not change much be- 
low the strange quark mass. Consequently, the chi- 
ral limit result is consistent with ao(1450). Further- 
more, it is emphasized that the would-be tj'tt ghost 
states give negative contributions to the ao correla- 
tor in both the quenched case [H[ and the partially 
quenched case 16[ when the quark mass is lower 
than the strange. Thus it is essential to take out 
these quenched or partially quenched artifacts before 
one can confidently obtain the ao mass. In the case 
of cr, a calculation of the tetraquark mesonium with 
the pseudoscalar-pseudoscalar 4-quark interpolation 
field [13] has been performed. From the deviation 
of the lowest state mass from that of the expected 
two-pion scattering state, it was suggested [20] that 
a bound mesonium state is seen. 

In the present work, we shall use the overlap 
fermion [211 ] to calculate ao and the scalar tetraquark 
mesonium. The overlap fermion has the benefit of 
having exact chiral symmetry at finite lattice spac- 
ing. Since chiral symmetry plays a pivotal role in 
these mesons in chiral effective theories (in partic- 
ular, it is concluded in the recent dispersion analy- 
sis of tttt scattering that the occurrence of a is on 
the basis of "current algebra, spontaneous symme- 
try breakdown, and unitarity" [22]]), we believe it is 
desirable to adopt a fermion action which explicitly 
exhibits the spontaneously broken chiral symmetry 
at finite lattice spacing. 

For the study of ao, it is well known that there 
are tj'tt ghost states in the quenched [H, EH and 
partially quenched [l6| cases. These ghost states 
become the physical t/7t and tj'tt states in the full dy- 
namical calculation without partial quenching. Such 
physical two meson states in full QCD have been 
seen when the quark masses are light enough [l5j . 
Thus, in order to obtain ao(1450) and possibly 
a (980) in the quenched approximation, one needs 
an algorithm which can fit multiple states includ- 
ing the tj'tt ghost states which lie lower in mass 
than the 1=1 scalar qq below the strange quark 
mass region. We have developed a sequential em- 
pirical Bay es method (SEBM) for constrained-curve 
fitting 23] to fit multiple states. This is based on 
the constrained-curve fitting of lattice data with 
Bayesian priors [24L |25l| . In the sequential empiri- 
cal Bayes method [231, one extracts the priors for 
the mass and spectral weight from a subset of data 
by fitting the two-point correlation function starting 
from the large time separation. First, one fits the 
ground state in a time window and then uses its fit- 



ted mass and spectral weight as priors to fit the first 
excited state in an extended time window. The pro- 
cess is repeated until time slices are exhausted. One 
then does a constrained-curve fit to the rest of the 
data set with the extracted priors. This method has 
been employed to fit the Roper resonance, the radi- 
ally excited nucleon, and (1535) on top of the rj'N 
ghost states [ll| . It has also been used in extracting 
radially-excited states of IP charmonium [261 ]. 

It turns out that, for the range of quark masses 
that we fitted, the 3-volume dependence (from a 
comparison of 16 3 x 28 and 12 3 x 28 lattices) of the 
spectral weights of the respective ghost tj'N two- 
particle scattering state and the bound one-particle 
baryon state come out in agreement with expecta- 
tion, as derived in Ref. [HIH^]. We regard this as a 
highly non-trivial test for the fitting method. Simi- 
larly, in the study of the pentaquark state O + (1540), 
we fitted the state in addition to the rj'KN ghost 
state and found from the volume dependence that 
it is in agreement with the KN scattering state 
for a large range of quark masses [13] ■ Through 
these studies, we are more confident that the fit- 
ting method is capable of fitting multiple states in- 
cluding the ghost states. It is of course limited by 
how good the data are and how many states one can 
fit, given the number of time slices of the lattice. 
We shall use this algorithm in the present study of 
ao(1450) and er(600). The smallest pion we have is 
182(8) MeV, which is substantially lower than most 
of the previous calculations of ao; this allows us to 
study the behavior of ao with the pseudoscalar me- 
son mass ranging from 1.3 GeV down to 182 MeV. 
This is important in revealing that the ao mass is 
very insensitive to quark mass in the range from 
strange quark down to the physical u/d mass. This 
turns out to have significant phenomenolo gica l im- 
plications on the pattern of scalar mesons [28| . For 
the study of the tttt mesonium with the tetraquark 
interpolation field on our lattice, it is crucial for the 
pion mass to be lower than ~ 300 MeV in order to be 
able to disentangle the tttt scattering states from the 
one-particle mesonium, in order to reveal the nature 
of the fitted states from the tetraquark correlator. 
Thus, for both the case of ao(1450 and c(600), it is 
essential to study them in the chiral regime where 
m-rr is smaller than 300 MeV. 



II. o (1450) AND 7fo(1430) MESONS 

Our calculation is based on data from 16 3 x 28 and 
12 3 x 28 lattices with 300 quenched Iwasaki gauge 
configurations [fi = 2.264) and overlap fermions 
with a lattice spacing a = 0.200(3) fm determined 
from f w (m w ). This makes our lattice sizes at 3.2 fm 
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and 2.4 fm, respectively. We will discuss the scale 
determination from the Sommer scale r$ later to as- 
sess the systematic error in scale setting. A subset 
of these quark propagators was used to study the 
quenched chiral logs in pion and nucleon masses |29l| , 
the Roper and 5*11 [ill ], and the pentaquarks [271 ] . 

There has been a concern that lattice spacing of 
0.2 fm may be too coarse for the overlap fermion 
and speculation that the range of the ove rlap Dirac 
operator may be as lon g as 4 lattice units [30| ■ How- 
ever, direct calculation [3l| at lattice spacings of 0.2, 
0.17, and 0.13 fm with Iwasaki gauge action reveals 
that the range of the operator is comfortably small in 
each of these cases (one lattice unit in Euclidean dis- 
tance and 2 units in "taxi driver" distance, the latter 



being defined as ttd = \ \ x ~ y\\i = YIl 



and it approaches zero toward the continuum limit. 
Thus, we don't think there is an issue regarding lo- 
cality of the overlap operator at 0.2 fm that we base 
our results on. 
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FIG. 2: ao correlators from the iptj) interpolation field for 
several quark masses with corresponding pion masses in 
MeV. 



We first report results on ao for which we use 
as the interpolation field. Since ao is an isovector, 
one only needs to calculate the correlator with the 
connected insertion. Shown in Fig. [5] are ao correla- 
tors as a function of time for 9 quark cases with pion 
mass from = 182(8) MeV to = 764(5) MeV. 
It is seen that for pion mass lower than ~ 600 MeV, 
the ao correlator starts to develop a negative tail, 



and it is progressively more negative at earlier time 
slices for smaller quark masses. This is a clear indi- 
cation that at least one of the ghost rj'n states, being 
lightest in mass, are dominating the correlator over 
the physical ao at larger time slices. This has been 
reported in the literature for the quenched [H, [l4| 
and partially quenched [lq ] calculations and the 
ghost contribution has been removed with the help of 
re-summed hairpin diagrams. The ghost rj'ir contri- 
bution in the ao correlator has been studied in the 
chiral perturbation theory [l3|, H(| . The one-loop 
hairpin diagram gives the following contribution to 
the isovector scalar meson propagator 



-Ar 2 m 2 



((p- k) 2 + ml)(k 2 +ml) 2 ' 



(1) 



where V3 = L x L y L z is the three-volume of the lat- 
tice, ro is the coupling between the scalar interpo- 
lation field and the rj and tt [T3 , IT(| or the matrix 
element {Q^tp^ 1 -k) , and = 2Nf Xt/f% is the hair- 
pin insertion mass which is related to the topological 
susceptibility Xt in the pure gauge theory. 

The Fourier transform (FT) of the dimensionless 
Ah(p) for the case of p = gives the following con- 
tribution to the scalar meson correlator on the Eu- 
clidean lattice 



G s (p = 0)=FT{a 2 A h (p = 0)} = - 
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where E n = y k 2 + m 2 and the (1 + E^t) factor 
is due to the double pole of the would-be rj' ghost 
propagator in the loop. N$ and Nt are the number 
of lattice points in the spatial and time direction re- 
spectively. The result for the more general partially 
quenched case has been derived in Ref. [HI and the 
corresponding expression for the would-be rf — N 
one-loop contribution to the nucleon correlator is de- 
rived in the study of Roper and Sii(1535) [TlT |. 
We shall use the expression 



W(l + E n t)e 



(3) 



where W is referred to as the spectral weight in later 
discussion, to model the fit of each ghost state con- 
tribution. Here we allow the energy, -Ey^, of the in- 
teracting r]'—w to be fitted to the data, but retain the 
double pole character of the prefactor 1 + E n t. This 
should be a good approximation when the rj' — n in- 
teraction is weak (N.B. in the large N c consideration, 
the meson-meson interaction is of the order 1/N C ) 
so that the prefactor due to the double pole in the 
would-be 7/ propagator remains largely valid when 
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higher orders are included. We have used a similar 
expression to fit the rj'N ghost states in the nucleon 
and £ii(1535) correlators [ll[ and found that their 
spectral weights have the correct 3-volume depen- 
dence for the range of quark masses we calculated. 
Thus, we shall employ this form to fit the clq corre- 
lator and, as a cross check, will examine if the l/E% 
dependence in Eq. ^ is borne out from the fit. 

We use the above-mentioned SEBM [23[ to per- 
form the curve-fitting with the weight W of Eq. ([3]) 
constrained to be negative, to reflect the ghost na- 
ture of the state as shown in Fig. [3J and the total 
energy of the would-be 77' and 7r constrained to be 
not far from the energy of the two non-interacting 
pions, i.e. 2 -y/Pn + m \ with discrete lattice momenta 
Pn = ^y 21 , n = 0, ±1, • ■ •. In the course of studying 
SEBM |23| , we extracted the priors from a subset of 
the data and applied them in a constrained fit of the 
rest of the data; we found that the results were very 
compatible, without detectable bias, with those ob- 
tained when the priors were applied to the full data 
set. In the present work, we extracted the priors 
from 100 gauge configurations and applied them to 
a fit of the remaining 200 configurations; we find that 
the results are very close to those obtained when the 
priors were applied to the full 300 configurations, ex- 
cept the latter had smaller errors. We shall report 
the results based on the 300 configurations here. For 
to^ > 250 MeV on the 12 3 x 28 lattice, we have been 
able to fit 4 states with the lowest and the third one 
being the ghost i]'tt states which are close to the non- 
interacting pair with each meson at zero momentum, 
and one unit of lattice momentum (pi), respectively. 
The second state has a positive weight and is inter- 
preted as the physical ao with the usual exponential 
form er Et . We show in Fig. [3] the fits to the ao 
correlators for several low quark masses. 

Due to the fact that there are expected to be two 
to three ghost states below the ao on the 16 3 x 28 lat- 
tice, we have, unfortunately, not been able to fit all 
of them to extract the physical ao with the SEBM fit- 
ting in this case. There are not enough time slices to 
fit 5 states. Otherwise, we could have compared the 
spectral weights from the 12 3 x 28 and 16 3 x 28 lat- 
tices and checked the expected volume dependence 
of the one-particle ao and two-particle 77' — 7r states. 
However, we could and did compare the ratio of the 
a correlators Ci2(t)/C\e(t) between the 12 3 x 28 
and the 16 3 x 28 lattices. As was derived in the 
course of studying pentaquark states [13] , the spec- 
tral weight of a one-particle state in the point-sink 
and point-sink correlator is proportional to unity; 
whereas, that of a weakly interacting two-particle 
state is proportional to I/V3. We show the ratio 
of the ao correlators Ci2(t)/Cie(t) in Fig. 0] for the 
cases of = 764(5) MeV and m„ = 182(8) MeV. 
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FIG. 3: Fitted ao correlators with a low-lying ghost rj'n 
state for several quark masses. 



In the case of the heavier pion, the ratio of the corre- 
lators for the whole time range is close to unity. This 
reflects the fact that there are no ghost states, so the 
lowest state is just the scalar qq meson and the ratio 
of the spectral weights is independent of the volume. 
On the other hand, the ratio for the = 182(8) 
MeV case is close to [V 3 (12)/V" 3 (16)]- 1 = 2.37 for 
t > 3, indicating that the lowest state is the ex- 
pected two-particle ghost rj'ir state. 
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FIG. 4: Ratio of ao correlators at 12 3 and 16 3 lattices for 
the cases of m n = 764(5) MeV and m m = 182(8) MeV. 
While the expected two-particle ghost t/'tt state shows 
volume dependence at lower quark mass, one particle ao 
state does not show any volume dependence at higher 
quark mass. 

In addition, we can check the dependence 
in the spectral weight W as suggested in Eq. ^j. 
We adopted the fitting form in Eq. © based on the 
premise that the higher loop diagrams are not im- 
portant. If the spectral weight indeed exhibits the 
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l/£^ dependence, it would lend support for such an 
assumption. We note that there are two ghost t/'it 
states below the physical ao on the 12 3 x 28 lattice 
for pion mass below 250 MeV. In this case, the low- 
est interacting would-be rj and 7r scattering state 
is close to the non-interacting pair with each meson 
at zero momentum; the second is close to the non- 
interacting pair with each meson having one unit of 
lattice momentum, i.e. p\. We have been able to fit 
4 states in the time window from t — 12 — 13 to 
t = 2. The third state has a positive weight and is 
interpreted as the physical ao with the usual expo- 
nential form e~ Et . Plotted in the left panel of Fig. [5] 
is the spectral weight W\ that we fitted for the lowest 
ghost rj — 7r state as a function of the pion mass. We 
see that as the m w decreases, it is quite singular. We 
fitted with 1 / E\ from the pion mass from 575 MeV 
down to 200 MeV. It is found that one can obtain a 
good fit down to ~ 270 MeV. Below ~ 270 MeV, 
there is a deviation. This is presumably due to the 
higher order effect. Similarly, we plot the spectral 
weight W2 in the right panel in Fig. [5l We see that 
it is non-zero below ~ 250 MeV and its magni- 
tude increases as E„ decreases. However, it does not 
increase nearly as fast as 1/E%. We conjecture that 
we may not have isolated the second ghost state in 
our SEBM fit and it may have included higher ghost 
contribution. Since the second ghost state starts to 
show up below = 250 MeV where the fitted spec- 
tral weight W\ starts to deviate from the expected 
1/m^ behavior, we do not quote results on ao below 
= 250 MeV. 
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FIG. 5: Fitted spectral weight W\ for the lowest ghost 
state (left panel) and W2 for the second ghost state in 
the ao correlator from the 12 3 x 28 lattice as a function 
of the pion mass and energy in GeV. 

Given the caveats of the ghost state fitting, we in- 
terpret the second state from our resultant fit from 
the 12 3 x 28 lattice which has the ordinary expo- 
nential form in time and positive spectral weight 
to be the physical ao- We plot its mass as a func- 
tion of the corresponding m\ in Fig. [6] together with 
that of a\ for comparison. The latter does not have 
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FIG. 6: ao and ai masses are plotted as a function of 
m \ . Also shown is the two pion mass (dashed line) which 
becomes lower than ao around the strange quark mass 
region. 



ghost state contamination and is thus easier to cal- 
culate. We see that above the strange quark mass, 
a\ lies higher than ao as expected from the quark 
model of heavy quarks. However, when the quark 
mass is smaller than that of the strange, ao levels 
off, in contrast to a\ and the other hadrons that 
have been calculated on the lattice. This confirms 
the trend that has been observed in earlier works at 
higher quark masses [HI, [TH, The chirally ex- 
trapolated mass ao = 1.42 ± 13 GeV suggests that 
the meson ao(1450) is a qq state. By virtue of the 
fact that we do not see ao(980) which lies lower 
than ao (1450), we estimate the spectral weight ra- 
tio (0|V^|a (980))7(0|y^Ml450)) 2 to be less than 
0.015 from the relative error of the correlator in the 
time window where ao is fitted. We also calculated 
the -fT (IT30) mass with the strange mass fixed at 
ma = 0.26833 which gives a vector meson mass cor- 
responding to the <j) mass and the u/d is extrapo- 
lated to the chiral limit. In this case, we also need 
to fit and remove the ghost Krj' states due to the 
hairpin diagram of the u/d quark. There is no ghost 
state due to the s quark, since according to Fig. [2J 
no ghost state is seen for quark mass greater than 
1/4 of that of the strange. As a result, we obtain 
the K (1430) mass at 1.41 ± 0.12 GeV and the cor- 
responding scalar ss state from the connected inser- 
tion to be 1.46 ± 0.05 GeV. Our findings are quite 
consistent with the experimental fact that if (1430) 
is basically degenerate with ao(1450) despite hav- 
ing one strange quark. This unusual behavior is not 
understood as far as we know and it serves as a chal- 
lenge to the existing hadronic models. 

It is known that the scale in the quenched ap- 
proximation is not determined uniquely. We note 
that if the Sommer scale r$ = 0.5 fm is used for the 
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scale, the lattice spacing will be 12% smaller, i.e. 
a = 0.175(3) fm. As a result, all the masses deter- 
mined above will be ~ 14% higher in the ro = 0.5 
fm scale. The same is true with the following calcu- 
lation of the 7T7T state and the tetraquark mesonium. 

Since we do not see ao(980) in the ao correla- 
tor with the tpjj} interpolation field, it leaves room 
for it to be something other than a qq state, e.g. a 
q 2 q 2 state as suggested in model studies. However, 
it is a challenge to verify it on the lattice due to 
the complication that there is a threshold KK state 
nearby (within 10 MeV). Therefore, we shall study 
the cr(600) first which, if present as a tetraquark 
mesonium, is several hundred MeV above the tttt 
threshold and several hundred MeV below the next 
7T7T scattering state with momentum close to pi . This 
is so, provided that the pion mass is lower than 
~ 250 MeV, and it may present the best hope of 
detecting such a state without the worry of entangle- 
ment with the collateral two-meson scattering states. 
Since the lowest pion mass in our case is 182 MeV, 
we are in a position to examine it. 



III. a(600) MESON 

The a meson was first postulated by M.H. John- 
son and E. Teller as a classical field to explain the 
saturation properties and binding energies of nuclei; 
they estimated a mass ~ 500 MeV from the sur- 
face energy (32|. It has been suggested that the a 
is partially responsible for the enhancement of the 
AI = 1/2 decay in K -> tttt [Hj]. Although it has 
been put back in the particle data table on account 
of the D + — > 7r~7r + 7r + experiment [34j . its experi- 
mental existence is still not fully settled due to the 
complication that its large width is as large as its 
mass. Recent dispersion analysis (35| using the Roy 
equation has produced a resonance pole in tttt scat- 
tering with high precision. The mass is given as 
441±g 6 MeV with a width of 544+^ MeV. Lattice 
QCD, in principle, is capable of resolving the issue 
about its existence. 

Resonance can be viewed as a mixture of a bound 
state and the scattering states in the usual poten- 
tial model description of scattering and resonance. 
In a coupled channel approach, one can couple a 
bound state in the continuum with the scattering 
states via a coupling potential resulting in a bound 
state leaking to the continuum with a shift in mass 
and acquiring a width. On a hypercubic lattice 
with periodic boundary conditions, the available mo- 
menta are p n = n^-, n = 0, ±1, ±2 ■ • • and therefore 
the one- and two-meson spectra for a certain quan- 
tum number are discrete. Imagine one considers a 
very large box where the two-meson spectrum has 



closely-spaced levels. When one looks at the spec- 
tral weights of the correlator in the channel with 
a specific quantum number, there will be envelopes 
of states with enhanced spectral weights which are 
the result of the mixture of a bound state and the 
nearby scattering states. They are the finite box 
representation of the resonances in the continuum. 
The "width" of the structure would reflect how far 
apart in energy the mixing takes place. If we start to 
decrease the size of the box, the spectrum is going to 
be less dense and states are further apart from each 
other. The number of states under the envelope di- 
minishes. When the box is small enough so that the 
scattering states are spaced enough apart that none 
is expected to lie under the envelope, then only the 
bound state remains and thus can be identified as 
such. By being "far enough apart" we mean that we 
can define the separation of the scattering states A 
to be several times greater than the width T of the 
structure, i.e. A ^> V. In this case, the bound state 
and the scattering states are not mixed and they 
will each have a different volume dependence in their 
spectral weights. Comparing the spectral weights of 
the same unmixed state in two lattice volumes is an 
effective way of revealing the one- or two-particle na- 
ture of the state [TJ H3 when the interpolation field 
projects to both the bound and scattering spectra 
in the correlator with a definite quantum number. 
This is our approach of identifying, through the vol- 
ume study of their spectral weights, both a bound 
tetraquark state (the a) and a two-pion scattering 
state which are reasonably well separated and un- 
mixed. 

Experimentally, a is a very broad state with a 
width of 544 MeV according to the recent dispersion 
analysis [35[ . It is natural to ask if it is ever possible 
to delineate it with a lattice calculation in Euclidean 
space. To answer this, we shall reverse the above dis- 
cussion on the existence of the bound state. Suppose 
one finds a bound state in addition to the scattering 
states which are outside the "width" of the state at a 
relatively small volume. To gain information about 
the scattering phase shift and hence the real width, 
we now increase the box size. As the box size is in- 
creased, the energies of the scattering states will be 
lowered since the value of the discrete momenta de- 
creases. When the scattering state above the bound 
state is lowered to within the range of the "width" , 
it mixes with the bound state (this actually defines 
the range of the width) and the two states avoid 
the level crossing. From the energy of the scattering 
state one can deduce the scatte ring phase shift us- 
ing the Liischer formula [36|, [33, 123, |39|, |4(| ■ This is 
valid for elastic scattering irrespective of how broad 
the state is. This is studied in detail in a 2-D lat- 
tice model and a spin model [39[ which illustrate 
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FIG. 7: The lowest three states from the scalar 
tetraquark correlator as a function of for m n from 
182 MeV to 250 MeV on the 12 3 x 28 lattice. The 
solid lines indicate the energies of the two lowest non- 
interacting pions in S-wave with lattice momenta po — 
and p± — 2-7T / 'La. 



how the scattering state mixes with the bound state 
and gives rise to the phase shift as the volume is in- 
creased. In a sense, by varying the lattice volume, 
hence the momentum, one can use a scattering state 
to mix with the bound state and scan the spectrum 
to obtain the phase shift and therefore the width of 
the resonance. The information of the width can 
also be obtained by determining how far apart in 
energy the scattering and bound state start to avoid 
the level crossing. 

In the present manuscript, we are only concerned 
about the existence of a and will leave the study 
of its width to the future. In this vein, we have 
chosen the lattice size such that the two lowest 7T7t 
scattering states are expected to be more than half 
of the experimental width, i.e. 272 MeV, away from 
the expected a mass at ~ 600 MeV. If this expected 
result is borne out, we can use the volume test to 
discern the particle content of the states and thereby 
distinguish the bound a from the two-pion scattering 
states. This is basically our strategy to seek the 
existence of a . 

To confirm the existence of cr(600) beyond doubt, 
one needs to identify both the tetraquark mesonium 
and the collateral irir scattering states. Secondly, 
one needs to work on a lattice where the scattering 
states and the bound state are well separated (e.g. 
further apart than half of the width of the "would- 
be" resonance) to avoid admixture; this is in order 
to discern the nature of these states separately to 
make sure that a is indeed a one-particle state, not 
a two-particle scattering state. To this end, we used 
the adaptive Bayes curve-fitting method |23j as de- 
scribed above to fit the ground state and the ex- 
cited states of the tetraquark correlator with a lo- 



cal interpolation operator ip"f5' t P' l l , "f5' , P f° r both the 
source and the sink. We note that it does not mat- 
ter what interpolation field one uses for the cal- 
culation, as long as it has overlap with the states 
of the corresponding quantum number. Being lo- 
cal, the pseudoscalar-pseudoscalar operator that we 
adopt will have vector- vector, scalar-scalar, • • • com- 
ponents after Fierz transform. This aspect has been 
discussed in the calculation of the exotic I — 2 pp 
states [4l| . It is also shown in the pentaquark study 
that different interpolation fields are related through 
Fierz transform [27j and the masses from these in- 
terpolation fields were verified to be the same in a 
lattice calculation [42| . Unlike in the ao correla- 
tor, there are no rfn ghost states to worry about 
in this tetraquark channel. As in Ref. [2(j, we con- 
sider only the connected insertion, not the single and 
double annihilation insertions. They are likely to 
preferentially project to the higher qq and glueball 
states 20]. To verify this, we use two point sources 
at t = and t — 8 and a zero-momentum wall 
source in the Coulomb gauge at t = 14 to calculate 
the disconnected-insertion correlator at time separa- 
tions t — to = 0, 6, and 14 and found that they are an 
order of magnitude smaller than the corresponding 
connected-insertion correlator. This shows that the 
annihilation diagrams are not likely to change the 
results of the connected insertion qualitatively. We 
present our results on the 12 3 x 28 lattice in Fig. [7] as 
a function of m\ for the pion mass range from 182 
MeV to 250 MeV. We have fitted three states so that 
we can trust the results of the two lower ones. The 
lowest one is about 100(20) to 60(10) MeV below 
the 7T7T threshold. This is most likely the interacting 
state of two pions with energy close to and below 
that of two non-interacting pions at rest, since the 
interaction is attractive in the I — channel. It is 
shown that the energy shift of two interacting par- 
ticles in a finite box can be related to the infinite 
volume scattering length below the inelastic thresh- 
old in a systematic 1 / L expansion 36] . In particular, 
for two spinless bosons with mass m at rest, one has 



AE = E - 2m = - 



47rao 
mL 3 



1 



(4) 



where c\ = —2.873 and C2 = 6.3752 from one-loop 
calculation [36[. However, it is pointed out [43| 
that there are would-be rj' hairpin diagram contribu- 
tions at the one-loop order for ttt: scattering which 
spoils the relation between AE and the scattering 
length ao in Eq. Q. It has L° and L 2 terms in 
the 1 = channel in addition to the leading 1/L 3 
term in Eq. Since our calculation is done in the 
quenched approximation, using the full QCD one- 
loop chiral perturbation formula [36( to extract the 
scattering length of 7T7t scattering from the energy 
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FIG. 8: The calculated energy shift of the lowest 1 = 
nn state AE as a function of m n L for the 12 3 x 28 lattice 
(left panel) and the 16 3 x 28 lattice (right panel). The 
lines labeled with 5 indicate the quenched one-loop chiral 
perturbation results for S = 0.10,0.15 and 0.18. The 
one-loop full QCD chiral perturbation prediction is also 
given for comparison. 



shift in the finite box channel is not applicable. We 
shall, instead, compare our results to that derived in 
quenched chiral perturbation theory. The quenched 
one-loop 7T7T scattering energy shift in the finite box 
which includes the hairpin diagrams has been de- 
rived 43]. The energy shift is 

AE = E - 2m w = AE tICC + AE onc - loop (5) 

where the tree-level result for the 7 = channel is 

-7 



AE tICC = 



4/2X3' 



(6) 



with / T = 132 MeV at the physical pion mass and 
the one-loop result is given 43j as 



AE' 



one— loop 



O 



(m„L) 3 



where 



m 2 /3 



16^ 2 /2 



(7) 



(8) 



Wc interpolate B^m^L) and Ao(m- K L) listed in 
Ref. [43| for the range of m^L appropriate for our 
calculation on the 12 3 x 28 and 16 3 x 28 lattices 
and plot AE for 6 = 0.10,0.15 and 0.18 which 
cover the range of S corresponding to the Witten- 
Veneziano formula for the rf mass and from the 
study of quenched chiral log in the pseudoscalar me- 
son masses 29| . The results are presented in Fig. [5] 
together with our data on the 12 3 x 28 lattice (left 
panel) and the 16 3 x 28 lattice (right panel) from the 
lowest state in our calculation which we believe is the 
two-pion scattering state. We see that our data are 



[V 3 (12)/V 3 (16)r 

[V 3 (12)/V 3 (16)] C 
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FIG. 9: Spectral weight ratio W12/W16 as a function of 
m% for the lowest state (filled circle in Fig. [7} and the 
next lowest state (filled square in Fig. . 



reasonably consistent with the one-loop quenched 
chiral perturbation calculation [43| for the range of 
8, i.e. from 0.10 to 0.18, especially for m^L > 2.8 for 
the 12 3 x 28 lattice and m^L > 3.5 for the 16 3 x 28 
lattice. This is true despite of the fact that we have 
not included the disconnected insertion in our cal- 
culation. This is consistent with our earlier finding 
that the disconnected correlator is about an order of 
magnitude smaller than the connected one at several 
time separations. 

We note that the tree-level energy shift in Eq. © 
is the same in quenched and full QCD. If the energy 
shift from the loop is small compared to the tree re- 
sult in both the quenched and full QCD, then the 
quenched energy shift in Eq. ([5|) would agree with 
the full QCD case in Eq. To this end, we plot the 
energy shift from Eq. with the scattering length 



,1=0 



7m^ 



A — from the tree-level in Fig. [8] and find 

that it is about a factor of two to three smaller than 
our data in the low m^L region. This shows that the 
loop contribution in the quenched case is enhanced 
compared to that in full QCD and not small com- 
pared to the tree part and thus one cannot apply the 
full QCD relation in Eq. ^ to obtain the scattering 
length. One can only compare the energy shift to 
that of the quenched chiral perturbation theory as 
was done in the last paragraph. 

The third state in Fig. [7| with a large error bar is 
about 1 GeV above the lowest state. The fact that 
it is higher than the energy of the non-interacting 
two-pion state (each pion with momentum p± = 524 
MeV), as indicated by the higher solid line, is an in- 
dication that the highest fitted state is always higher 
than the true state as it inevitably includes the un- 
fitted higher states and hence cannot be taken as a 
good signal for a definite state. 

One interesting aspect of the spectrum is that 
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there is an extra state between the lowest irir scat- 
tering state and the third state which presumably 
encompasses the higher scattering states. It has a 
sizable spectral weight, as large as that of the lower 
7T7T scattering state. The mass is around 600 MeV. 
It is tantalizing to identify it with <r(600). To verify 
this, we study the volume dependence of the spec- 
tral weight of these states. It was advocated in the 
study of the Roper resonance, the pentaquark, and 
the ghost state [ll], |27| that one efficient way of dis- 
tinguishing a one-particle state from a two-particle 
scattering state in a finite box is to study the vol- 
ume dependence of its spectral weight. From the 
normalization factor of l/^/V^ for a particle in a box 
and the way the correlator is constructed, i.e. with 
a point source and a zero-momentum sink, the spec- 
tral weight of a one-particle state does not explicitly 
depend on volume; whereas, the spectral weight of a 
weakly interactin g tw o-particle state has an explicit 
I/V3 dependence [23]. This is true when the one par- 
ticle state is reasonably separated from the scatter- 
ing states so that the mixing of the two states is not 
strong enough to spoil the characteristic volume de- 
pendence of their respective spectral weights. Since 
we have two lattices with sizes 12 3 x 28 and 16 3 x 28, 
the spectral weight ratio for a two-particle state 
should be W 12 /W 16 = V 3 (16)/V 3 (12) = 16 3 /12 3 = 
2.37. Plotted in Fig. [5] are the ratios of the spectral 
weights for the lowest state in Fig. [7| and the first ex- 
cited state around 600 MeV. We see that the spectral 
weight ratio W\ijW\s for the lowest state clusters 
around 2.37, confirming our speculation that it is the 
interacting two pion state. On the other hand, the 
spectral weight ratio W12/W16 of the excited state 
near 600 MeV turns out to be consistent with unity. 
This suggests that this state is a one-particle state, 
not a 7T7T scattering state, and strongly supports the 
identification of it to be the <r(600). Furthermore, by 
virtue of the fact that the ratios of spectral weights of 
these two low-lying states are consistent with unity 
and V3(16)/V3(12), it suggests that the mixing of the 
two states, if any, is small. In other words, the two 
states with an energy separation of ~ 300 MeV are 
reasonably well separated compared to half of the 
decay width of <r(600). Extrapolated to the chiral 
limit, the mass of the one-particle state is 540 ± 170 
MeV. 

There have been concerns that the zero mode con- 
tributions which are finite volume effects may con- 
taminate the results when the volume is small for 
the quark mass under study. In the case of the 
chiral condensate, the zero modes can be avoided 
by doubling the contribution from the chiral sector 
which does not have zero modes 44]. For the pion 
mass calculation, the zero modes can be removed 
by considering the correlator of the pseudoscalar- 



pseudoscalar (PP) and scalar-scalar (SS) combina- 
tion (PP) - (SS) [Hill. We have studied the zero 
mode contributions to the pion mass calculation on 
the same set of lattices used in the present study [29[ 
by comparing the pion mass from the (PP) correla- 
tor and the (PP) — (SS) correlator and found that 
there is no difference in the pion masses within er- 
rors. We concluded that the volume is large enough 
for this lattice that the zero mode contribution is 
less than the error for the range of quark mass stud- 
ied. In the present study of a (1450), ifg(1430) and 
cr(600), we don't see any indication of mass diver- 
gence for small quark masses. We thus believe that 
the zero mode contributions are within errors in the 
connected insertions. 

Although lattice study of tetraquark states started 
some time ago [2(J lU, |4j H§| , we believe this is the 
first time that both the one-particle state and its 
concomitant scattering state are identified and their 
nature verified through the volume dependence of 
the spectral weights. Short of such identification, 
we don't think one is able to rigorously confirm the 
existence of the tetraquark mesonium. 

Finally, we note that the calculation of tetraquark 
state with the ip"/^M'Ys'>p interpolation field has first 
been attempted [20]. It was found that in the 1 = 
channel, the ground state is lower than the expected 
two interacting pion state from chiral perturbation 
theory and is interpreted as a bound state — a 
tetraquark mesonium. However, in this analysis the 
full QCD formula in Eq. ^ was employed which 
we pointed out earlier is not applicable to quenched 
calculations. If the quenched chiral perturbation cal- 
culation [43[ , which has a different lattice length de- 
pendence, is used one might come to a different con- 
clusion. 



IV. CONCLUSION 

To conclude, we calculated the isovector ag with 
the t/jip interpolation field and er(600) with the 
tetraquark interpolation field. With the overlap 
fermion, we have come down in the chiral region with 
very low pion mass (182(8) MeV) in the quenched 
approximation. After removing the fitted rfir ghost 
states, we found the lowest ao at 1.42±0.13 GeV and 
Kq at 1.41 ±0.12 GeV which are consistent with the 
experimental ao(1450) and Xq(1430) being the qq 
states and confirms the earlier findings in quenched 
and partially quenched calculations at higher quark 
masses. In addition, we have been able to fit the 
1 = scalar tetraquark correlator and have iden- 
tified, through the volume study of their spectral 
weights, both the lowest interacting two-pion state 
and a one-particle state at 540 ± 170 MeV. This sug- 
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gests that <r(600) does exist as a particle and it is 
a tetraquark mesonium. This is consistent with the 
recent dispersion relation analysis of tttt and KK 
scattering with the Roy equation which has led to 
the cr pole at 4411^ MeV with a width L CT = 544±|| 
MeV [35| . Further lattice calculations with dynami- 
cal fermion in the chiral region with < 300 MeV 
are needed to check these results. 
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